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Abstract 

A new class of completely integrable models is constructed. These models are defor¬ 
mations of the famous integrable and exactly solvable Gaudin models. In contrast with 
the latter, they are quasi-exactly solvable, i.e. admit the algebraic Bethe ansatz solu¬ 
tion only for some limited parts of the spectrum. An underlying algebra responsible for 
both the phenomena of complete integrability and quasi-exact solvability is constructed. 
We call it ’’quasi-Gaudin algebra” and demonstrate that it is a special non-Lie-algebraic 
deformation of the ordinary Gaudin algebra. 
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1 Introduction 


A quantum model is usually called exactly solvable if all solutions of its spectral problem 
can be found algebraically. Recently a new important class of the so-called quasi-exactly 
solvable models has been discovered [||, [|], ||, These models are distinguished by the fact 
that only certain limited parts of their spectra admit an algebraic construction. A detailed 
exposition of different theories explaining the phenomenon of quasi-exact solvability and 
proposing constructive methods for building and solving such models can be found in papers 
0 §■ 0 ! 0 0 and in the book & 

The quasi-exactly solvable models can be considered as deformations of exactly solvable 
ones [|I!|]. Usually, the deformation of a given exactly solvable model leads to an infinite 
sequence of quasi-exactly solvable models with different hamiltonians and different number 


of exact solutions |12]. 

For example, the simplest exactly solvable model of a multi-dimensional harmonic oscil¬ 
lator with hamiltonian 




i —1 


dxl 


2 2 


+ 0 b~x 
i— 1 


( 1 . 1 ) 


can be deformed into an infinite sequence of quasi-exactly solvable models of multi-dimensional 
sextic anharmonic oscillators with hamiltonians 


H n = ~ E ^2 + E + 2ar2 


i= 1 


i=l 


d 


J2 b i x i - 2n - 1 - - 


-i =1 


+ a 2 r 6 , n = 0,1,. 


( 1 . 2 ) 


In contrast with Ql.lj ), each of models (1.2) has only (n+d)\/(n\dl) algebraically constructable 
solutions j|]. Here the role of a deformation parameter is played by a. 

It is worth stressing that both models (1.1) and (|L^) are completely integrable in the 
sense that they have enough number of mutually commuting integrals of motion (see e.g. 
|l2| ). Thus, the above example clearly demonstrates that the deformation procedure may 
preserve the integrability property. There are many other examples of such a sort which can 
be found in ref. |T|]. 

It is natural to assume that any completely integrable and exactly solvable quantum 
model can be deformed into an infinite set of integrable and quasi-exactly solvable models. 
If this assumption is true, then the following interesting problem immediately arises. 

The Problem. It is well known that most of physically interesting integrable and exactly 
solvable models can be obtained and solved in the framework of the celebrated R-matrix 
approach [13, 14, ^]. A mathematical structure lying in the ground of this approach and 
responsible for both the phenomena of integrability and exact solvability is the so-called Yang 
Baxter (YB) algebra (see e.g. the book jl(|). Using its generators, one can easily construct 
not only all the integrals of motion of a model (and thus ensure its integrability) but also 
build all eigenvectors and eigenvalues of the corresponding spectral problem (and thus ensure 
its exact solvability in the framework of the so-called Algebraic Bethe Ansatz Ju|). Let 
us now consider any quasi-exactly solvable deformation of such a model that preserves its 
integrability property. If such a deformation does exist, then it is natural to ask ourselves 
if there is any algebra responsible for both the phenomena of integrability and quasi-exact 
solvability of the deformed model? In other words, is there any analogue of the YB algebra 
(some kind of ”quasi-YB algebra”) whose generators could be used not only for constructing 
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the integrals of motion of the arising quasi-exactly solvable models but also for building their 
algebraically constructable eigenvalues and eigenvectors? We mean here some analogue of 
the Algebraic Bethe Ansatz. And the last and, may be, the most important question: is it 
possible to interpret this quasi-YB algebra as a deformation of the ordinary YB algebra? 


It is obvious that if solution of this problem does actually exist and the quasi-YB alge¬ 
bras will be found, then, developing a regular method for building their representations we 
obtain a big factory of new integrable and quasi-exactly solvable models. Since the ordinary 
YB algebras give rise to many physically interesting exactly solvable systems of quantum 
mechanics, statistical physics and field theory, it would be natural to expect that the set 
of quasi-exactly solvable models associated with quasi-YB algebras also will contain a lot of 
physically interesting systems. 


It seems that the posed problem has a positive solution. In this paper we consider 
the simplest case of the so-called sZ(2)-Gaudin algebra G[sl( 2)] which is closely related to 
the classical YB algebra and leads to the completely integrable Gaudin models solvable by 
means of Bethe ansatz (section |2j) . In sections ||| and Q] we construct a special non- 

Lie algebraic deformation of this algebra which automatically produces an infinite series of 
integrable models. We call this algebra the ”quasi-Cr[sZ(2)] algebra” and the corresponding 
models - the ”quasi-Gaudin models”. In sections 11 and [7] we show that all these models 
are quasi-exactly solvable in the sense that they admit only a partial Bethe ansatz solution of 
spectral problem. In section [|| we construct a simple realization of quasi-G[sZ(2)] algebra and 
later, in section |9j, demonstrate in an independent way that quasi-Gaudin models associated 
with this concrete realization are actually quasi-exactly solvable. In the last section |10 


we 


introduce a concept of quasi-sZ(2) algebra which can be viewed as a special limiting case of 
quasi-G[s/(2)]. 


2 The Gaudin algebra 


In this section we remind the reader some basic facts concerning Gaudin algebras, their 
representations and properties of the associated Gaudin models. More detailed exposition of 
this subject can be found in refs. 

The Gaudin algebra G[sZ(2)] is an infinite-dimensional extension of the ordinary sl( 2) al¬ 
gebra. Its three generators, S( A) = |5 _ (A), S' 0 (A), S'" 1 " (A)} with A G C obey the commutation 
relations 


17, 18 


5°(A)5°(/i) - S°(ji)S°( A) 
S ± (X)S ± ( f i)-S ± ( f i)S ± (X) 

S°(X)S ± (n) - S ± (fi)S°(X) 
S-(X)S + ^)-S + ^)S~(X) 


0 , 

0 , 

_ L 5 ± (A)-5 ± (/i) 
H — A 

o 5°(A )-5°(/i) 
H- X 


generalizing those of sl( 2). Using (2.1), one can easily prove that the operators 


C(A) = S°(X)S°(X) - ±S~(X)S + (X) - ±S + (X)S~(X) 


( 2 . 1 ) 

( 2 . 2 ) 
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form a commutative family, 


[C(A),C(/i)]=0, (2.3) 

and thus, can be interpreted as integrals of motion of a certain quantum completely integrable 
model. The latter is known under name of Gaudin model. 

The role of the “Hilbert space” in which the operators C'(A) act is played by the repre¬ 
sentation space of Gaudin algebra. In order to construct it one needs to fix the lowest weight 
vector |0) and the lowest weight function /(A) obeying the relations 

5°(A)|0)=/(A)|0), S-(A)|0)=0. (2.4) 

After this, we can define the representation space as a linear hull of vectors 

I£i,...,£ m ) = S + (£ m )...S + (6)|0) (2.5) 

with arbitrary m and £i,...,£ m . We shall denote this representation space by 
The “Schrodinger equation” for the Gaudin model reads now 


C(A)0 = </>€W fW . (2.6) 

The beauty of this equation lies in the fact that all its solutions can be obtained algebraically. 


For example, using commutation relations (2.1) and formulas (|2.4| ) it is easy to check that 
the lowest weight vector is always a solution of the Gaudin problem: 


C(A)|0) = (>(A) + ^/(A)) |0). 


(2.7) 


The remaining solutions of equation (2.6) can be obtained by means of the so-called algebraic 
Bethe ansatz 


d> = s + (U)s + (Zm- 1 ) • • • 5 + ( 6 ) 5 + ( 6 )| 0 ), 


( 2 . 8 ) 


in which m is an arbitrary non-negative integer and £i,... ,£ m are some unknown numbers. 
It can be demonstrated that vector (2.8) is an eigenvector of the Gaudin operator with the 
eigenvalue 


E(x) = /2(A)+ Jd (A)+ £ /<A) m 


i=l 


A -£i 


if the following conditions hold 


E 

k=l 


1 


, ti ~ ik 


+ /(&) = 0, i = 1,... ,m. 


(2.9) 


( 2 . 10 ) 


These conditions are known as Bethe ansatz equations. It is known that the constructed 
solutions with m = 0,1,... expire all possible solutions of the Gaudin spectral problem and 
therefore the latter is exactly solvable. 
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3 The quasi-Gaudin algebra 


In this section we consider a special modification of Gaudin algebra which also leads to 
completely integrable quantum systems. 

Let S n (A) = {5-(A),5°(A),5+(A)}, n £ Z, A £ C denote the operators obeying the 
relations 


S° n (X)S^) - S^S^X) = 0, 


0/,,\ e0/ 


5± ±1 (A)5±(^)-5,f ±1 (/r)5,f(A) 

5° ±1 (A)5±(/r)-5±(^)^(A) 


= 0 , 


S~ +1 (\)S+ (n) - S'+_ 1 (/i)5“(A) = 2 


S±( A) - S±(n) 
fi-X 

s° n (\)-sM 


jJL — A 


(3.1) 


We consider (3.1) as the defining relations of a certain infinite-dimensional algebra. It is not 
difficult to see that this algebra is very similar to the Gaudin algebra but, in contrast with 
the latter, it is not a Lie algebrafj We call it the “quasi-Gaudin algebra” or, more concretely, 
“quasi-G[s/(2)]”. 


The algebra (3T) has many remarkable properties. First of all, only from the “quasi¬ 
commutation relations” ([0]) it immediately follows that the operator valued functions 


C n ( A) = ^(A)5°(A) - ^S- +1 (A)S+(A) - ^+_ 1 (A)5“(A) 


form commutative families for any n: 


[C'n(A), C n (n)\ = 0. 


(3.2) 


(3.3) 


These functions are obvious generalizations of the Gaudin integrals of motion. Note, however, 
that the commutators between C n ( A) and C m { A) do not vanish if m ^ n. 


The property (343) suggests to interpret C n { A) as generating functions of commuting 
integrals of motion for certain completely integrable quantum models. In fact, ( |3.3| ) describes 
an infinite sequence of such models, since the integer n can be chosen arbitrarily. This is the 
first important difference with the Gaudin case when an analogous construction leads to a 
single completely integrable model — the Gaudin model. 

Exactly as in the case of Gaudin algebra, the construction of these models needs the 
specification of the “Hilbert space”. It would be natural to identify it with the representation 
space of our algebra. The latter can be constructed in the same way as in the Gaudin case. 
First of all, we need the lowest weight vector |0) which, as in the Gaudin case, should be an 
eigenvector of both the lowering and neutral operators S~( A) and 5°(A). The most general 
relations expressing this fact and compatible with the quasi-commutation relations ([0]) read 


S°(A)|0> = (F(A) +nG(A)) |0), S"(A)|0> = nG(A)|0>, 


(3.4) 


2 In section |§] we demonstrate that (jhl]) can be viewed as a deformation of the Gaudin algebra. 
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where F( A) and G( A) are certain arbitrarily fixed functions^. Note that the lowest weight 
vector |0) is not annihilated by the “lowering operators” S~(X) except for the case of n = 0. 
This fact, which will play a determining role in our further considerations, gives us the second 
important difference with the Gaudin case. 

After fixing the “lowest weight functions” F( A) and G( A), we can define the representation 
space as a linear hull of vectors 


l£i,...,{ m >*=S+ ({ m )S+ —l+fc(6n—1 )...s+ (? 2 )s+ ( 6 )| 0 ) 


(3.5) 


with arbitrary k, m and £i,..., £ m . The convenience of this dehnition comes from the fact 
that, due to the relations ( 0 , the vectors ( |3.5| ) are symmetric with respect to all permuta¬ 
tions of numbers £i,..., (exactly as in the Gaudin case). We denote the representation 


space of algebra (0) by W>(A),G(A)- 


4 The quasi-Gaudin model 


Now we are ready to construct the integrable models associated with algebra (3.1). 
postulate that Schrodinger equations for these models read 


We 


Cn{X)(f) n — E n (X)(j> n , (j) n £ Wp(X),Q(X ), n — 0,1,.... (4.1) 

So, formula (|4.lD dehnes an inhnite sequence of integrable models, each of which is completely 
determined by a triple {F(X), G(X), n}. 


5 A trivial solution 

It is known that in the Gaudin case the lowest weight vector is always a solution of the 
Gaudin spectral problem. Is this true in the case of models (4.1)? In order to answer this 
question, we should simply act by the operator C n {X) on the vector |0) and look at the result. 
For this it is convenient to rewrite the operator C n {X) in a little bit different form. We use 
the formula 


d 


S~ +1 (X)St(X) - 5+_i(A)5-(A) = -2—5°(A), 


(5.1) 


which follows from the quasi-commutation relations (3T) in the limit n —» A and, substituting 
it into the expression (|3.2[) for C n ( A), obtain 


(5.2) 


Cn( A) = S0(A)S°(A) + ^S° n ( A) - S+_i(X)S~ (A). 


From (5.2) and formulas (pAh it immediately follows that 


C'n(A)lO) = {(F(A) + nG(A )) 2 + ^ (F(A) + nG(A))] |0) + nG(A)5+_ 1 (A)|0). (5.3) 

3 Stictly speaking, the G( A) functions in the first and second formulas of (^0|) may differ from each other 
by a certain n-dependent constant factor. However, using the transformations 5? (A) — > 5° (A), £“( A) — > 
CnS~ (A), S+ (A) — > c n+i'S'il'fA) which do not change the commutation relations ([hll), we can simply remove 
this difference. 
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Now it becomes clear that vector |0) could be a solution of problem (4T) only if n = 0. In 
this case we obtain the relation 


C 0 (A)|0) = (f 2 (A) + ^F(A)) 0) 
which is similar to the Gaudin relation (2.7). 


(5.4) 


6 The Bethe ansatz 


In analogy with the Gaudin case, let us try to find the solutions of spectral equations (4.1) 
by means of the algebraic Bethe ansatz. 

It is natural to try to take the Bethe vector in the form 


<t>n = S+ +k (£ m )S+ -1+k (fm-l) • • • Sj +t ({2)Sj +1 Ki)|0) 
with some k and m, and, using the relations 

S^lMSn-lW - S+-lWS%-l(ri 


Cn(X)S:_M - A) = 2- 


A — // 


( 6 . 1 ) 


( 6 . 2 ) 


which follow from the basic relations (3.1), try to transfer the operator C n (A) to the right. 


From formula (G2) it is seen that, in order to start performing such a permutation, it is 
necessary to take 


m + k = n — 1. (6.3) 

Now note that each permutation of the operator C n (A) with raising generators decreases its 
index by one. This means that after m permutations when this operator appears in front 
of the lowest weight vector, it will have the index n — m. The standard presciptions to 
Bethe ansatz technique imply that the hamiltonian of an integrable system, appearing after 
all permutations in front of the lowest weight vector, should be absorbed by it. But such 
an absorbtion is possible only if the lowest weight vector is an eigenvector of a hamiltonian. 
According to the result of previous section, the lowest weight vector is an eigenvector of the 
operator C n _ m (A) only if 


n — m = 0. 


(6.4) 


Comparing formulas ( |6.3D and 6G we can conclude that the only case when the ansatz (HJ 
for equation ( 0 ) may lead to some algebraic solutions corresponds to the choice 


k = —1, m = n. 


(6.5) 


7 The Bethe ansatz solution 

After using the restrictions ( |6.5D , the ansatz (|0]) takes the form 

K = S^ n )S^_ 2 (^_r)... 5+(e 2 )5 0 +(e 1 )|0). (7.1) 

Let us now check that it actually contains solutions of the problem ( 0 ). This can be 
demonstrated exactly in the same way as in the ordinary Gaudin case. 
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Repeating the reasonongs of ref. |17|, we transfer the operator C n { A) to the right by 


using the quasi-commutation relations (3.2). The neutral generators of quasi-Gaudin algebra 
appearing after such a transference also should be transfered to the right. Exactly as in the 
case of a (7-operator, any permutation of a S'°-operator with raising generators forming the 
Bethe vector decreases its index by one. Finally, after completing all possible transferences, 
we obtain a lot of terms with operators Co (A), Sq (A) and ,Sq (G) standing in front of the 
lowest weight vector |0). After this one should get rid of these operators by using formulas 
(3.4) and (|5.4|). The result has the form 


C n (A)S'+_ 1 (^ n )... S 0 +(£i)|0) = A(A,6, • • •, £n)SU(£n) • • • S+foM + 

E ... Ete+i)E-i(A)E- 2 te-i) • • • E( 6 )|o), 


(7.2) 


i=l 


where A(A,£i,...,£ n ) and Bi(£ i,..., £ n ), i = 1,..., n are some explicitly constructable func¬ 
tions. We see that the right hand side of (7.2) consists of two parts. The first part is 
proportional to the Bethe vector (j> n , while the second one consists of n terms which do not 
have such a form (these are the so-called ” unvanted terms”). This means that the coefficient 
function A(A,£i,..., £ n ) determine the eigenvalues of the operator C n (X) provided that the 
unwanted terms are absent. But in order to get rid of the unwanted terms it is sufficient to 
require the vanishing of their coefficients -Bj(£i,..., £ n ) = 0, i = 1,..., n This leads us to a 
system of n equations for n unknowns £i,..., £ n , which is nothing else than a typical system 
of Bethe ansatz equations. The explicit form of these equations reads 


n 


i 


k=i,k& ^ £ k 


+ F(€i) — 0, i — 1, 


,n, 


(7.3) 


and the final expression for the eigenvalues E n ( A) is given by the formula 

E n (X) = F\ A) + ^F(A) + E F{X) x ~_ F ^ i) 


(7.4) 


i=l 


Formulas 0, 0) and (fF4|) complete the solution of problem ( flip. I t is not difficult to 
see that this is only a partial solution since the linear hull of vectors (7T) with fixed n gives 
us only a certain negligibly small part of the whole representation space defined in section 
|H The above consideration enables one to conclude that the models which we obtained are 
typical quasi-exactly solvable models. 


8 Simplest realizations of quasi-Gaudin algebra 

In previous sections we did not discuss the realizations of algebra ( |3.ip . Now it is a time to 
do this. In fact, we know only one realization which can be constructed from the generators 
of standard Gaudin algebra S( A) satisfying the commutation relations 0- In order to 
construct it, we will need the special limiting operators 


S ± = lim A5' ± (A), 

A— KX) 
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5° = lim AS 0 (A), 


( 8 . 1 ) 













which form the sl( 2) algebra. Using these operators, let us take 

n + / - 5° 

A — a 

n + f-S° + b 
A — a 

n + f-S° + 2b 
A — a 


Sn (A) = 5“(A) + 
S°n( A) = 5°(A) + 
5+(A) = 5+(A) + 


( 8 . 2 ) 


where a and b are complex parameters and / = lim.\_ >0 o A/(A). It can be easily checked by 
direct calculations that operators actually satisfy the quasi-conmrutation relations (3.1). 


The lowest weight functions F( A) and G( A) characterizing the representation of this al¬ 
gebra read 


F( A) = /(A) + 


A — a 


G( A) = 


1 


A — o 


(8.3) 


Now it is absolutely clear that algebra (|0]) is a deformation of the Gaudin algebra. The 
role of the deformation parameter is played by number a. If a —* oo, then the n-dependence 
of the operators (|8.2| ) disappears and they transform into ordinary generators of Gaudin 
algebra. Respectively, the quasi-commutators in ( |3.1|) become the ordinary ones. As to the 
formulas (|8.3| ) defining the representations of algebra o, they, in the limit a —» oo also 


transform into the relations (2.4) for Gaudin algebra, 


9 Quasi-exact solvability: a different point of view 


Substituting flS.2|) into (|3.2|) , we obtain the form of the operators 

CnW = S°(A)5°(A)-i5-(A)5 + (A)-i5 + (A)5-(A) + 

2S°(X)(n + b + f — S°) — S~(X)(n + 2b + f — S°) — 5+(A)(n + / — S°) 


A — a 


+ 


b(b-l) 

(A - a)2 


(9.1) 


Let us now check in an independent way that the operators (0) actually describe quasi- 
exactly solvable models. First of all, note that if the formulas ( |S.2| ) hold, then the repre¬ 
sentation spaces of Gaudin and our algebras coincide. Denote by the linear hull of all 
vectors S + (£k)... S' + (£i)|0) with arbitrary £i,... ,£*. and k < n. It is known that if /(A) is a 
rational function, then dim<L n < oo, for any n. From the obvious relations 5 ± (A)<J> ri C 
5°(A)<I> n C <!>„, (n + / — S , °) ( I> n C <L n _i and (n + / — S°)<I> m C <l> m for m / n it immedi¬ 
ately follows that the operator C n ( A) admits only one algebraically constructable invariant 
subspace, 4> n . This subspace is finite-dimensional and therefore the models (4.1) are quasi- 
exactly solvable. 


10 Conclusion. From quasi-G[s£(2)] to quasi-s/(2) 

It is known that the sl( 2) algebra can be considered as a limiting case of the s/(2)-Gaudin 
algebra s The generators of the latter are given by the formula ( |8.l[) . 


9 












It is natural to ask what kind of algebra will appear if we consider an analogous limit of 
the deformed algebra ( [3.1D ? In analogy with the Gaudin case, we can define the generators 
of this algebra by formulas 


St = lim AS±(A), 

A—xx) 


5° = lim XS° n (X). 

A—»oc 


( 10 . 1 ) 


Multiplying the quasi-commutation relations (|3.l|) by A/r and tending both A and /r to infinity 
we easily derive the relations between these generators, which read 


qO q0 q0 q 0 _ rj 


SnilSn - St ±1 St = 0, 


c0 c0 _i_ 

*^n=bl*^n ^n » 


Sn+lS+ - St^S~ = 25°. 


( 10 . 2 ) 


We can consider ( 10. 2|) as dehnig relations of a certain modification of the sl( 2) algebra which 
can be called “quasi-s/(2) algebra”. Despite the fact that it is not a Lie algebra, it has many 
properties similar to those of the ordinary sl(2). For example, it has a quasi-analogue of the 
Casimir operator, 


_ q 0 q 0 1 Q— Q+ ^ Q+ Q- 

'-'n o^Ti+l^n r> l^n 


2 2 

which quasi-commutes with all generators: 

/~< c+ _ q + o— _ q — ri q 0 _ q 0 

— D n— ]Dn— 1; ^n^n- |_l — ‘- 5 n -|_iC-'n+l, O n O n — J n ^ n . 


(10.3) 


(10.4) 


The representations of quasi-s/(2) algebra can be constructed in the same way as in quasi- 
Gaudin case. Defining the lowest weight vector |0) by the formulas 


S°\0) = (F + nG)\0), S~\0) = nG\0), 


(10.5) 


where F and G are certain arbitrarily fixed numbers, we can define the representation space 
Wf,g as linear hull of vectors 


I m)k = St +k St_ 1+k ...S+ +k S+ +k 


| 0 >, 


( 10 . 6 ) 


with arbitrary k and m. If we are interested in the spectrum of quasi-Casimir operator, then 
it is easy to see that, in contrast with the standard sl( 2) case, it is not infinitely-degenerate 
and contains only a few number of exactly constructable eigenvectors. For example, the 
only possibility for vector ( 10.6 ) to be an eigenvector of Q is realized when k = — 1 and 
m = n. In this case, the corresponding eigenvalue is equal to F(F — 1). So, one can say that 
the quasi-Casimir operators represent the simplest (and, in some sense, trivial) quasi-exactly 
solvable models. 

It is remarkable that the generators of quasi-sZ(2) algebra can be realized in terms of 
the generators of ordinary sl( 2) algebra. The corresponding formulas can be obtained after 
substitution of formulas (|__ ^_ 

(10.7) 


]) into (|8.1|) and read 
S~ = S~+n + f-S°, S° n = n + b + f, S+ = S + + n + 2b + f - S°. 


This realization corresponds to the choice F = f + b and G = 1 in (|10.5 ). We do not know 
at the moment if there are other, more general, realizations of quasi-,s/(2) algebra. 
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Concluding this paper, let us mention two problems which arise as immediate conse¬ 
quences of the results obtained above. These are: 1) construction of general quasi-Gaudin 
algebras and models and 2) construction of general quasi-Lie algebras. However, one should 
stress again that our main goal is to develop a theory of general quasi-YB algebras. The 
first step in this direction will be done in our next publication which will appear soon in 
hep-th archive and in which we intend to construct quasi-exactly solvable deformations of 
the inhomogeneous XXX spin chain and investigate the associated quasi-Yangian Y[^(2)]. 
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